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We study charmonium properties below and above T c up to 1.8T C , on quenched anisotropic lattices. Information 
of the spectral functions is extracted using the maximum entropy method and the constrained curve fitting. We 
also calculate the color singlet and averaged free energies and evaluate the charmonium spectrum with the potential 
model analysis. The relation between the lattice result of the spectral function analysis and the potential model 
is discussed. 



1. Introduction 

To investigate the properties of quark gluon 
plasma (QGP) in heavy ion collision experiments, 
theoretical prospects are significant, since such 
processes include complicated interactions among 
large number of particles. Changes of charmo- 
nium states have been regarded as one of the 
most important probes of plasma formation pQ, 
because the potential model calculations predict 
the mass shift of charmonium near T c 2., and 
J/ip suppression above T c [3]. However, lattice 
QCD simulations have indicated that the ther- 
mal properties of hadronic correlators are much 
more involved than weakly interacting almost free 
quarks ^ . Recent studies of spectral functions of 
charmonium suggest that a hadronic excitation 
of c-c system may survives above T c |5I6I7| . This 
seems to conflict with the predictions of potential 
model approaches and a naive picture of QGP. 

In this study we discuss the above disagree- 
ment between the spectral function analysis and 
the potential model. On one hand, we extract 
the information of the spectral function from the 
temporal charmonium correlator. On the other 
hand, we also extract the static quark potentials 
from the color singlet Polyakov loop correlation 
and Wilson loop. The result of potential model 
calculation using the latter is compared with the 
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former. 

These calculation are performed on quenched 
anisotropic lattices. This enables simulations at 
temperatures from T ~ (Nt = 160) to 1.75T C 
(Nt — 16) keeping the lattice cutoff unique, to 
avoid uncertainties caused by the cutoff depen- 
dence in, for example, self-energy contribution 
to the potentials. The anisotropic lattice is also 
useful to keep the number of data points of the 
temporal correlators sufficiently large, which is 
important for detailed analysis of the spectral 
functions. We extend the simulation reported in 
Ref. 5 which was performed at the renormalizcd 
anisotropy a a ja T = 4, the spatial lattice cutoff is 
a^ 1 = 2.0GcV, and with the plaquette gauge and 
O(a) improved Wilson quark actions. The critical 
temperature almost corresponds to N t = 28. 



2. Spectral functions 

To extract the spectral function from the tem- 
poral correlators, we adopt the maximum entropy 
method (MEM) and the constrained curve fitting 
(CCF) The result of the former provides the 
prior knowledge required for the latter. Combin- 
ing the two methods, the results may be more re- 
liable and quantitative than with one of them. In 
our previous study |5| we found that MEM fails to 
extract the spectral function from the correlator 
of local operators at T > 0. Therefore we adopt 
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spatially extended operators to enhance the low 
frequency mode of the spectral function. How- 
ever the smeared operators may lead to an artifi- 
cial peak, and thus careful analysis is necessary to 
distinguish the physical results from the artifact 
ones. 

Although our MEM result shows that the spec- 
tral functions have peak structure in PS and V 
channels (corresponding to r\ c and J/ip) at all 
temperature, the difference of result for differ- 
ent smearing functions, we call "smeared" and 
"half-smeared" , exists at higher temperature, es- 
pecially above 1AT C . This means that the peak 
structure of the spectral function at high tem- 
perature might be artificial. Furthermore we find 
large default model function dependence of the re- 
sults, in which the position of the peak is stable 
but the peak width has large dependence. There- 
fore it is difficult to study the spectral function 
quantitatively using only MEM. 

CCF analysis is performed based on the result 
of MEM, assuming the form of fitting function as 
the Breit-Wigner type, 
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where Ri, m,i and 7, are overlap, mass and width 
of the i-th peak respectively, and N ter m — 2 is 
used in this analysis. The result of lowest peak 
position (mass) and width are shown in Fig^ 
Unfortunately, our results of the CCF is unstable 
against the prior knowledge as inputs. Therefore 
the systematic uncertainties will be larger than 
the quoted statistical error in Fig|T] 

In spite of large systematic uncertainty, we can 
find the following tendency. Below T c , the spec- 
tral function has a peak with almost the same 
mass at T = and vanishing width. Above T c , 
the peak stays at similar position as at T < 
while its width tends to grow as temperature 
increases. As pointed out in the MEM analy- 
sis, however, the difference between the results of 
smeared and half-smeared correlators shows that 
a disappearance of the peak structure at higher T 
cannot be excluded, especially in the V channel. 
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Figure 1. Temperature dependence of lowest peak 
mass, TO;, and width, 71, from the result of CCF 
analysis. 

3. Static quark free energies 

In this section we compare the results of spec- 
tral function with the potential model calculation. 
Here we define the following potentials. 

V ave (r) = -ln(TrL(r)TrLt(0))/(iV t a T )(2) 
V sing (r) = -ln(TrL(r)tf(0))/(N t a T ) (3) 
Vwuson(r) cx -\nW(f,t)/t (4) 

where L(r), W(f,t) are Polyakov and Wilson 
loops, respectively. V ave is adopted in the for- 
mer potential model calculations. V S i ng can not 
define in gauge independent form, then we calcu- 
late it in the Coulomb gauge, in which the V s i ng 
will be equivalent with a gauge invariant defini- 
tion |5] . We can also define the other color singlet 
potential using the Wilson loops, V w u son , which 
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Figure 2. Color singlet static quark potentials 
y s in 9 and Vwtison, below and above T c . T = 
result is also shown in both figures. 

is calculated with the same way as at zero tem- 
perature. V S ingiet is thermal average of the energy 
spectrum of the color singlet q — q system, while 
Vwiisonioop is the lowest one of the spectrum. 

Figure [21 shows the result of V S i ng , Vwason be- 
low and above T c respectively. Above T c , unfortu- 
nately, we fail to extract the Vwason and present 
only Vsing in the bottom panel. We note that our 
results of V aV e are consistent with previous stud- 
ies, and draws the same conclusions as Refs. |2I3| . 

We solve the Schrodinger equation using the 
quark potential parametrized to the form, V(r) = 
Ai + A2/r+A-3-r+A4-r 2 . Table^Jsummarizes the 
result for the differences of binding energies from 
T = 0. Below T c , the CCF analysis found the 
mass difference of —(0 — 10) MeV. Corresponding 
results of the potential model with V s i ng devi- 



Table 1 

Result of the Potential model analysis. The num- 
bers mean difference from T = result. 
T/T c {N t ) V smg (MeV) Vwason (MeV) 

0.88(iV 4 = 32) +(71 - 87) -(0 - 2) 
0.97(iV t = 29) +(50 - 58) -(21 - 23) 
1.08(iVt = 26) -(170 - 180) none 
1.40(iV t = 20) unbound none 



ate from the CCF results, while those of Vwason 
agree. Above T c , V S i ng can hold a bound state up 
to 1.08T C , while at higher temperatures cannot, 
or at least the wave function become too broad 
compared to the spatial lattice size. 

In summary, the disagreement between the re- 
sults of spectral function analysis and the poten- 
tial model is qualitatively absent if one uses the 
color singlet free energy as the static quark poten- 
tial. This is the same conclusion as the Bielefeld 
group ^U] . Detailed analysis of the quark poten- 
tial indicated that the potential from the Wilson 
loop most well reproduces the result of spectral 
function analysis below T c . 
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